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Abstract 

Large liquid argon TPC detector programs such as LBNE and LAGUNA-LBNO will be able to make measurements of the proton 
lifetime which will outperform Cherenkov detectors in the proton decay channel p —> K^v. At the large depths which are proposed 
for such experiments, a non-negligible source of isolated charged kaons may be produced in the showers of cosmogenic muons. 
We present an estimate of the cosmogenic muon background to proton decay in the p —> K^v channel. The simulation of muon 
transport to a depth of 4 km w.e, is performed in the MUSIC framework and the propagation of muons and secondary particles 
through to a cylindrical 20 kt LAr target is performed using Geant4. An exposure time of 100 years is considered, with a rate of 
< 0.0012 events/kt/year at 90% CL predicted from our simulations. 
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1. Introduction 

An attractive framework in which to embed the Standard 
Model (SM) of particle physics is a Grand Unified Theory 
(GUT) in which each of the three independent gauge coupling 
constants of the SM symmetry group SU(3)c ® SU(2)l ® U(1)y 
are unified at a high-energy scale, Agut- 

Any GUT model of hadrons, leptons and gauge interactions 
would necessarily imply the violation of baryon-number con¬ 
servation at the Acux-scale m. If one pursues such a the¬ 
ory, the fundamental observables of Acux-scale physics would 
be closely related to the stability of the proton. In particular, 
the lifetime of the proton and the dominant proton decay prod¬ 
ucts could indicate a preference towards a specific GUT model 
and could also provide an insight into physics below the Acux- 
scale. 

In the simplest extension to the SM, one can extrapolate 
the three coupling constants to high energies such that the dif¬ 
ferent couplings become the same order of magnitude above 
10'^ GeV 121. One caveat with this procedure is that the cou¬ 
plings do not unify at a single energy scale 011 as one might 
expect from the most simple GUT 0. The observed stability 
of protons in the decay mode p —> 0 presents another 

problem for the simplest GUT model. These problems can be 
avoided in extended GUT models by suppressing or forbidding 
this decay mechanism. Extended GUT models can incorporate 
supersymmetry 0 in which the favoured process for proton de¬ 
cay can be p —» K^v 0, whilst allowing the SM couplings to 
unify at a single energy 0. 
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Amongst the next generation of high precision proton-decay 
detectors will be large liquid argon (LAr) time projection cham¬ 
bers (TPC) such as LBNE or LAGUNA-LBNO 0[IO). These 
experiments will introduce LAr fiducial volumes of the order of 
tens of kt at depths of at least 1.5 km. LAr TPC imaging will 
have a higher efficiency to measure the mode p ^ K^v com¬ 
pared to Cherenkov detectors, and as such they will outperform 
existing detectors such as Super Kamiokande IfTTl . Given the 
high efficiency for particle identification in LAr TPC detec¬ 
tors, the largest background to the p K^v signature will be 
from cosmogenic muons and neutrinos. 

In this paper, we present the results of our simulation of the 
number of backgrounds to this signal from mesons pro¬ 
duced in the showers of cosmogenic muons. We compare our 
results to the study performed by Bueno et al in Ref. IfT^ . Com¬ 
pared to the latter study, we perform the first full Monte Carlo 
simulation of particle production, transport and detection that 
includes cosmogenic muons and all secondary particle^ We 
also consider a depth and scale of the detector volume which is 
more applicable to the LBNE and LAGUNA-LBNO detectors. 

2. Simulation framework 

In our model, a LAr TPC detector is positioned at a depth 
of 4 km w.e., which corresponds to the proposed depths of the 
LBNE and LAGUNA-LBNO detectors. We model a total mass 
of 20 kt of LAr, which is close to design specification of the 
LAGUNA-LBNO detector. The detector is modelled as a cylin¬ 
der of LAr with a diameter of 30 m and a height of 20 m. The 


*We define ‘primary’ particles as those which are present at the surface at 

the Earth, and ‘secondary’ particles as those which are produced at any depth 
below this. 
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Figure 1: The distribution of muon energies on the top and side surfaces 
of the cuboid described in the main text, as sampled and simulated by MU- 
SIC/MUSUN [T4]|T3. a total of lO’ muons are considered in this distribution. 

LAr cylinder in encased within a stainless steel container with 
a thickness of 5 cm. One can compare this to the study per¬ 
formed by Bueno et al IfT^ which reports a prediction of the 
background rate using a 100 kt volume of LAr at a depth of 

3 km w.e. simulated in the FLUKA lfT3l framework. 

The designs of large-scale LAr detectors include up to 2 m 
of non-instrumented LAr which will separate TPCs from the 
physical walls. This region of non-instrumented LAr defines 
a ‘wall’ for this study; as no interaction or energy loss occur¬ 
ring in this volume can be seen by the TPC and is unlikely to 
be caught by a light detector, unless a very sophisticated light 
detection system is used. We do not model this region of the 
detector, and instead treat this boundary as one of the physical 
walls of the detector. 

In this paper, the simulation of particle propagation is per¬ 
formed in two stages. In the first stage, only muon transport is 
considered and the interactions of secondary particles are ne¬ 
glected. In the second stage, all particles including secondary 
particles are fully simulated. 

In the first stage of the simulation, muons are propagated 
from the surface of the Earth through a vertical depth of 

4 km w.e. using the MUSIC transport code nma. A simple 
model of standard rock with proton number Z = 11, nucleon 
number A = 22 and density p = 2.65 g cm“^ is used in the 
muon transport simulation. The overall muon flux normalisa¬ 
tion is provided by the Gaisser parameterisation m , which has 
been modihed for large zenith angles. The muons are then sam¬ 
pled on the surface of a cuboid using the MUSUN code oa. 
The cuboid has a height of 32 m and horizontal dimensions of 
40 m X 40 m, within which the detector cylinder is centred. 
The energy spectra of a sample of 10^ muons which have been 
fully transported by MUSIC/MUSUN to the top and sides of 
the cuboid surface are shown in Figure [T] 

The second stage of the simulation is performed using 
GEANT4.9.(f]n7l. in which all primary and secondary particles 
are transported from the surface of the cuboid until all surviving 


^ The Geant4.9.6 ‘Shielding’ physics list is used, and the muon-nuclear 
interaction process is additionally switched on. 


particles have exited the cuboid volume. In the second stage of 
the simulation, all particle information and associated energy 
depositions in the LAr volume are recorded. 

3. Event selection 

In this section, we present event selection criteria which 
maintain high efficiency for selecting p ^ K^v signal events 
whilst rejecting cosmogenic muon-induced backgrounds. 

Kaons from free proton decay are expected to have a kinetic 
energy of 106 MeV. After accounting for the Fermi motion of a 
proton in a nucleus, the nuclear binding energy, and subsequent 
re-scattering of the kaon, the kinetic energy of a kaon emitted 
from a nucleus of argon is expected to be smeared across the 
range 0 to 200 MeV lITSl . In a LAr TPC detector, the kinetic 
energy measurement of charged kaons is further broadened by 
the energy resolution, which we estimate following the proce¬ 
dure presented in Ref. lfT9l . In order to retain > 99% of the sig¬ 
nal events we require that the energy deposition from a charged 
kaon due to ionisation and scattering, < 250 MeV. 

By constraining the total energy of all proton decay products 
to the proton mass, we require that sum of and the energy 
deposition of subsequent decay products, should be less 
than 1 GeV. Meanwhile energy deposition from other particles 
not associated with the kaon, its interactions or decay products, 
Eother, in the same event should be smaller than 50 MeV; other¬ 
wise this additional energy deposition would be clearly visible. 

We reject charged kaon tracks that originate from within 
10 cm of a detector wall to ensure that the candidate signal 
events from proton decay are not mismatched with the events 
originating outside of the detector. This results in a decrease 
in efficiency of 2%. We also require that all tracks are not 
within 10 cm of the detector walls, to ensure for proper energy 
evaluation and particle identification. Events containing such 
partially-contained tracks are rejected. The requirement for no 
activity within 10 cm of the detector wall is shown to reject 87% 
of events with no primary muons in the hducial volume. 

At the simulated depth of 4 km w.e., we hnd that muons pass 
through the LAr detector at a rate of 0.078 s ', so the average 
time between the two muons crossing the detector is approx¬ 
imately 12.8 s. Assuming that the maximum drift time, and 
hence the duration of the event record, will be about 10 ms, the 
probability of a muon crossing the detector within any 10 ms 
time window will be approximately 10“^. Rejecting events 
that contain a muon, where the muon’s track length is greater 
than 20 cm for clear identification, results in the reduction of 
p —> K^v detection efficiency of about 0.1%. 

LAr TPC detectors can identify charged kaons in the range 
of interest (< 250 MeV) via their dE/dv in LAr lfT9ll . Eor low 
energy charged kaons the identification efficiency is assumed to 
be 100%, which will lead to a slightly conservative estimate of 
the background to p —> K^v. In order to fully account for the 
kaon identification efficiency, detailed studies will need to be 
performed for each specific detector design. Although it is as¬ 
sumed that this method can identify charged kaons with a high 
efficiency, using dE/dx provides no separation between and 
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K states. The following procedure will provide some degree 
of separation. By considering the dominant /T* decay modes, 
one would expect to find a in 95% of decays, either di¬ 
rectly from the or via a subsequent decay. In these decay 
modes, and K can be separated in the event of either or 
capture in the liquid argon. Before considering the muon 
capture lifetime, the procedure of identifying positive kaons as 
having a decay chain which features yu* —> e* would reject 95% 
of negative kaons. After comparing the muon capture lifetime 
in argon lIMIl to the muon lifetime ED, the K rejection factor 
is reduced to 82%. In addition to decay chains featuring 
yU* ^ e*, we also permit events featuring (where X 

is any other set of particles) in order to account for the remain¬ 
ing 5% of decays. 

To summarise, we propose the following selection criteria for 
p ^ K^v events in large LAr TPC detectors; 

1. There is exactly one in the event. 

2. There are no muons with track length > 20 cm in the de¬ 
tector volume. 

3. There is no activity within 10 cm of the detector wall. 

4. The decay chain includes the decay e* or 

for electrons with Eg > 5 MeV. 

5. The total energy deposited by the excluding decay 
products, satishes < 250 MeV. 

6. The total energy deposited by the and by decay 
products, satishes (E^^^ + E^'^) < 1 GeV. 

7. The total energy deposited in the rest of the hducial vol¬ 
ume satishes Tsothei- < 50 MeV. 

4. Results 

In total, an exposure time of 100 years has been simulated. 
The number of events that are expected to pass the sequen¬ 
tial selection criteria described in Section]^ is presented in Ta- 
ble[2 After applying all cuts, we do not observe a single muon- 
induced event in our simulations. Figure shows the dis¬ 
tribution of E^^'- and Fother for events with charged kaons after 
applying cuts 1 - 4, as specihed in Section The region re¬ 
tained after the respective selection criteria for both and 
Fother are applied, is bound by the blue-dashed lines. It is found 
that even if the requirement on the kaon kinetic energy is loos¬ 
ened such that E^^'- < 400 MeV, there will be no events in the 
region of interest. 

We hnd that the main source of events events is due to 
the production of outside of the LAr TPC volume. For 
events passing the selection, the ionisation energy peaks be¬ 
low 200 MeV, but these events are rejected by requiring that the 
total energy deposition not associated with the charged kaon is 
lower than 50 MeV. 

The absence of a single muon-induced event from our 
event selection is converted to an 90% Confidence Level (CL) 
upper limit of 0.024 background events per year in a 20 kt liquid 
argon detector, or 0.0012 events/kt/year. This result is consis¬ 
tent with earlier estimates HI. The presented event selection 
reduces the efficiency for p K^v selection by less than 2%, 
which is dominated by the requirement for no activity within 
10 cm of the detector wall. 
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Figure 2: A scatter plot of the total energy deposited by mesons by ionisa¬ 
tion versus the total energy deposition not associated with the kaon in the event. 
Kaons are selected by applying cuts 1 - 4 as specified in Section]^ The region 
retained after the respective selection criteria for the two quantities is bound by 
the blue-dashed lines. There are no events observed with Mother < 50 MeV. The 
data represents the statistics after 100 years of exposure. 


5. Conclusions 

We have presented an estimate of the muon-induced back¬ 
ground to proton decay in LAr TPC detectors. We consider the 
background to the process p K^v, which is favourable in 
extended GUTs models. 

In our simulation, a cylindrical LAr TPC detector is posi¬ 
tioned at a depth of 4 km w.e., which corresponds to the depths 
of the proposed LBNE and LAGUNA-LBNO detectors. We 
consider a total mass of 20 kt of LAr, which is close to design 
specification of the LAGUNA-LBNO detector. 

We hnd that the main source of muon-induced background 
events is due to muon-induced mesons produced outside 
of hducial volume. We show that background events can be 
rejected by constraining to the kinematics of the proton decay 
as well as requiring that the total energy deposition in the event 
not associated with the is less than 50 MeV and that there is 
no activity within 10 cm of the detector wall. 

After considering an exposure time of 100 years, we set a 
90% CL upper limit on the number of background events as 
0.0012 events/kt/year. Furthermore, the presented event selec¬ 
tion is estimated to reduce the signal efficiency by less than 2%, 
which is dominated by the requirement for no activity within 
10 cm of the detector wall. 
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Table 1: The number of expected events after 100 years of exposure, as a function of sequential selection criteria and the K- production mechanism. Events are 
accepted into the table if the total energy deposition in the event is less than 2 GeV. Numbers presented in brackets indicate the subtotal number of K~ events. 
N refers to nucleons in the interaction. 


K* parent 

Exactly one K- 

No muon 

No activity near 
wall 

(P* OtK*) -7 

e* +X 

< 250 MeV 

< 

1 GeV 

£'other < 50 MeV 

fi* +N 

255 (43) 

59(10) 

0 

- 

- 

- 

- 

71* +N 

134 (20) 

79 (14) 

3(0) 

3 

3 

3 

0 

(pip) + N 

13(2) 

7(1) 

0 

- 

- 

- 

- 

y + N 

8 

6 

0 

- 

- 

- 

- 

S- 

I 

1 

0 

- 

- 

- 

- 


118 (28) 

63 (15) 

31(2) 

31(0) 

24 

24 

0 

{njn) + N 

11(1) 

9(1) 

0 

- 

- 

- 

- 


12(2) 

10(1) 

0 

- 

- 

- 

- 

Total 

552 (96) 

234 (42) 

34 (2) 

34 (0) 

27 

27 

0 
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